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ABSTRACT 

A photon or neutrino gas-semi-contained by a baryonic species through scattering- 
comprises a rather pecuUar MHD fluid where the magnetic field is truly frozen only 
to the co-moving volume associated with the mass density. Although radiative dif- 
fusion precludes an adiabatic treatment of compressive perturbations, we show that 
^ . the energy equation may be cast in "quasi-adiabatic" form for exponentially grow- 

' ing, non-propagating wave modes. Defining a generalized quasi-adiabatic index leads 

, to a relatively straightforward dispersion relation for non-axisymmetric magnetoro- 

' tational modes in the horizontal regime when an accretion the disk has comparable 

00 ' stress contributions from diffusive and non-diffusive particle species. This analysis is 

, generally applicable to optically thick, neutrino-cooled disks since the pressure contri- 

(N . 

butions from photons, pairs and neutrinos, all have the same temperature dependence 
whereas only the neutrino component has radiative heat conduction properties on the 
C~| . time and length scales of the instability. We discuss the energy deposition process and 

C^h' the temporal and spatial properties of the ensuing turbulent disk structure on the 

Q \ basis of the derived dispersion relation. 

H , 
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1 PRELIMINARIES 

Understanding the magnetorotational instability is key to the development of realistic models of accretion. As a physical 
process, the MRI justifies two sine qua non ingredients of accretion disk theory: entropy generation from the differential 
shear flow and turbulent angular momentum transport. Consequently, clarifying the dynamics that leads to MRI initiated 
turbulence in a few relevant astrophysical regimes is paramount in uncovering how the accretion process takes place. 

In the context of accretion onto compact objects, the stress associated with a radiative particle species often dominates 
the dynamics at high temperatures but our understanding of its effect on MHD processes is incomplete at best. Agol & Krolik 
(1998) pointed out that compressive radiation-MHD wave modes will be strongly suppressed by diffusive loss of pressure 
support in a range of wave numbers: Cph/c < < 1, where k = imip k represents the wave number normalized to the mean- 
free-path and Cph is the phase speed of compressive waves in the single fluid view. Such analysis is applicable to accretion disks 
at large wave numbers where the effect of inertial forces and stratification are minor. For the MRI, however, both inertial and 
stratification effects play key roles in dictating the dynamics and polarization of the modes. 

Blaes & Socrates (2001) have reported on a radiation-MHD dispersion relation for axisymmetric modes that comprises the 
MRI under rather general circumstances. Our analysis is different from theirs mainly in the consideration of a purely toroidal 
field where the compressibility of the non-axisymmetric modes leads to the largest impact from radiative diffusive damping. 
Furthermore, aiming to understand ultra-relativistic flows very near black holes, we emphasize toroidal mode dynamics 
because these modes may become predominant when shear effects begin to overwhelm other inertial accelerations leading to 
fastest growing modes at large azimuthal scales (Araya-Gochez 2002). The scale and degree of intermittency in the turbulent 
flow bear a direct impact on gravitational wave emission models from hyper-accreting black holes (Araya-Gochez 2003) . For a 
relativistically hot accretion disk with comparable vertical pressure support from diffusive and non-diffusive species (see §2), we 
solve for an approximate dispersion relation involving Lagrangian displacements in the horizontal regime {^e/^r — kr/ke 0). 
This simpliflcation affords a lowering of the rank in the dispersion relation to second order (in uj'^) which, in turn, enables 
a relatively straightforward approximate treatment of the effects of radiative heat conduction (in contrast with Blaes and 
Socrates' more complex but exact, relation for axisymmetric modes). 
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2 RADIATIVE VISCOSITY VS DIFFUSION 

The half optical depth to scattering of a standard radiation-pressure-dominated a-disk is related to the vertical scale height, 
Ti, and to the angular rotation profile, n(r), through 

Tdisk = (1) 

where Oort's A parameter is normalized A = ^dinrlnfi (= —3/4 for a non-relativistic, Keplerian flow). This relation ensues 
from balancing the local heating rate, oc MO^A, with the local cooling rate from radiative diffusion, oc Vrp; while eliminating 
the accretion rate in favor of the vertically integrated component of the viscous stress responsible for angular momentum 
transport oc a Hp. Eq 1 is thus explicitly sensitive only to the local nature of the cooling (although it is implicitly subject to 
a suitable vortical gradient of heat deposition; see, e.g., Krolik 1999). 

In the neutron-rich, neutrino cooled gas associated with a hyper-accreting black hole (see, e.g. Popham et al. 1999), 
the pressure contributions from radiation, pairs and neutrinos, all have the same temperature dependence, oc a/6 T*, with 
relative contributions varying only by internal degrees of freedom times particle statistics factors: 2x1,4x7/8, and 6 x 
7/8, respectively (with only one helicity state per neutrino). Thus, the neutrino pressure is never greater than ~ pthormai 
where pthcrmai ~ ll/12ar* includes the pressure from photons and pairs. Yet, only the former species has radiative diffusion 
properties on a dynamical time scale. Thus, in this setting the neutrinos are the lone diffusive species and their stress is 
identified with prad below (e.g., radiative = diffusive species throughout this paper). 

Partial pressure support by the non-diffusive components (baryons, toroidal field and the thermal component) does not 
modify Eq 1 as long as the contribution to the disk's flare from this extra pressure is inversely proportional to the mass density 
decrease from the same, e.g., ID expansion. Consequently, for a neutrino cooled a-disk, a replacement in the source of opacity 
by neutrino scattering in the non-advective accretion regime: 0.1 M© sec~^ < M < 1 M© sec~^ (Di Matteo et al. 2002), is all 
that is needed for Eq 1 to hold. 

On the other hand, the "parallel" size the of fastest growing MRI eddies is /eddy = ^MFii — ^au/^- ^'^^ magnetic angular 
momentum transport such that a ~ (^'Ait/Cw)^' the isotropic diffusion time through these eddies, tdis ~ , turns out to be 
similar to the timescale for fastest MRI modes to develop tmni — —-^otrt (^®® §5)- Moreover, the optical depth through these 
eddies is (with = ptot / p) 

Teddy ~ X Tdisk ~ j — — » 1- (2) 



c. 



n 



2Av^ 



Since Compton/neutrino drag go 6(3 and k ~ 0{ac^/c^), these estimates indicate that the effect of 

photon/neutrino viscosity is negligible for the optically thick eddies of a radiatively cooled disk. 



3 A LAGRANGIAN FORMULATION OF THE MRI IN SOFT MEDIA 

We take a semi-local approach by computing instantaneous (co-moving) growth rates in the limit kr/kg and referring to 
these as the fastest growing modes on the implicit understanding of their transient nature (since kr{t) = fcor — [dinr^^] k^t). 
In terms of a local, co-moving observer, azimuthal wave numbers are no longer discrete and consequently, neither are the 
CO- moving frequencies (Ogilvic & Pringle 1996). Further, we discard field curvature and radial gradient terms. A simple 
meridional stratification profile sets the physical scale-length of the problem: d^ hip = l/H, with gas, radiative and magnetic 
pressures tracking the unperturbed density profile. 

The first order relation between Lagrangian, A, and Eulerian, S, variations (see, e.g., Chandrasekhar & Lebovitz 1964, 
Lynden-Bell & Ostriker 1967) is 

A = S + $-V. (3) 
The Euler velocity perturbation, v = 5V, is thus related to the Lagrangian displacement, ^, through 

V = {9t V • V} e - (^ • V)V ^ ia^ - e^.rl^- (4) 

The algebraic relation follows from the assumption of differential rotation and from writing exp i{u;t + mip + kzz) dependencies 
for ^, while a = lo + mfi denotes the co-moving frequency of the perturbations. 

These geometrical equations have their more traditional equivalents in the so-called shearing sheet approximation where 
a co-moving, "locally Cartesian frame" (r,Lp,6) {x,y,z), is used along with the linearized shear velocity field, V(x) = 
[dinrf^] a;lj,, to treat the problem locally while introducing the coriolis terms by hand. (The reader is invited to examine a 
companion paper: Araya-Gochez (2002), where the inertial terms are induced in a fully covariant, geometric way). 

The linearized equations of motion for the displacement vector correspond to the Hill equations (Chandrasekhar 1961, 
Balbus & Hawley 1992) 

-2Qi^ = -4Aoortn + ^ , i^ + 2Qir = ^, and G = ^ , (5) 
where each over-dot stands for a Lagrangian time derivative and f, for the sum of body forces. 
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The Lagrangian perturbation of mass density and the Eulerian perturbation of the field follow straightforwardly from Eq 
[3] , and from mass and magnetic flux conservations (under ideal MHD conditions, notwithstanding the peculiar nature of the 
fluid under study) 

^ = -V-C, <5B = Vx(^xB). (6) 
Using the latter relation, the Lorentz force in the co-moving frame is readily laid out 
''i^'^''^ = C X { V(V.O+V|^-VBV(l6.C)-lgVB(V.O} 

^-^■^ -<u X {{ki^i - feB^B)k + kU - l^k^h^i-} (7) 

where 1^ is a unit vector in the direction of the unperturbed field and where the scalar operator Vg = 1^ • V. Note that the 
term {ki^i — fcgfg) = k±^± may be interpreted as a restoring force due to the compression of field lines (Foglizzo & Tagger 
1995). 

Next, we look in detail at the non-magnetic stress terms. 

In the single fiuid view, the Lagrangian variation of the specific pressure gradient, /0~^Vpr+g, contains two terms (see, 
e.g., Lynden-Bell & Ostriker 1967): one oc Ap~^, and another one oc AVpr+g. In terms of the displacement vector, the first 
term is proportional to the equilibrium value of Vpr+g which is negligible in the local treatment ((x radial gradient when 
^e/^r —^ $)■ For the same reason, the Eulerian and Lagrangian variations of the specific pressure gradient are nearly identical. 

The "thermodynamic" pressure term is thus given by 

ivp.+g) =r^V(V.O ^-^^ c^kfc,^,, (8) 

where, for a heterogeneous fluid, T{= d[inp]lnp) represents a generalized adiabatic index (see, e.g., Chandrasekhar 1939, 
Mihalas & Mihalas 1984). 

Assembling f from Eqs [7 & 8] , and plugging this form into Eqs [5] yields 

^ + 2 n X ^ = {c? k + *;^„(k - fcg lg)}(k • ^) + vl^ikl i - fee^B k) - AAQ^r Ir- (9) 

This equation agrees with the matrix dc-composition of Foglizzo & Tagger (1995) for F = 1. Moreover, these authors also 
pointed out that because of vertical stratification and rotation, the polarization of the slow MHD (e.g., MRI) modes obey 

when kr/ke — > 0. Eqs [9 & 10] impose an anisotropy constraint on the components of the Lagrangian displacement 

''^^^ "~fT2e ''ii^i' ^-^'^ii^ih (11) 

where e = p^/p,+g. 

In this horizontal regime, with kz(,z finite, the dispersion relation derived from Eq [9] and from the above constraint reads 



where all frequencies are normalized to the local rotation rate, = 4(1 + A) is the squared of the epicyclic frequency, and 
= (k- VAif )/n is a frequency related to the component of the wave vector along the field (in velocity units). This expression 
matches the form given by T. Foglizzo (1995). 

In terms of the function, A(r, 0), we find the fastest growing wave numbers to be given by 

g|=-2A + (^)x|-^|, where C s 1 + (i-A)A + ^1 + (1 - A)A, (13) 

and where the expression in the curly brackets corresponds to the growth rate of the modes. 

Notably, the compressibility of non-axisymmetric, MRI modes is imprint on deviations of A from unity 

^ = (1-A)(ifc||)^||, (14) 

e.g., the degree of compression gets stronger with field strength and, naturally, with a softer equation of state. 



4 A QUASI-ADIABATIC INDEX FOR WAVE PHENOMENA 

In a heterogeneous fiuid composed of a neutral baryonic plasma (with ideal gas properties) plus a neutrino and/or photon 
component, ideal MHD dictates that the field lines remain frozen to the charged massive component while pressure perturba- 
tions associated with such material gas must track density perturbations on all scales. On the other hand, coUisionaJ coupling 
between the gas and the diffusive species translates into a scale dependent containment of the pressure perturbation associated 
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with the radiative species. Consequently, acoustic wave modes partially supported by radiative stress will be damped from a 
non-adiabatic loss of pressure support when the diffusive and wave time scales are comparable (Agol & Krolik 1998). Such 
phenomenology can be quantitatively implemented by considering the mathematical equivalent to the equation of state for 
the diffusive component of the fluid. 

Wc start with the co-moving frame, frequency-integrated, radiative transfer equation-corrected to first order in the fluid's 
motion, (3, and with isotropic, elastic scattering as the lone source of opacity: 

dtl{n) + n ■ V/(n) = (1 + 3/3 ■ n) J - 2/3 ■ H - (1 - n • /3)J(n). (15) 

A tilde means normalization to the mean free path, i^tp, or to the mean crossing time, imtp/c, while /i/(n) is the radiative 
intensity in the direction of n, and J and H are the first two moments of the intensity. 

Agol and Krolik (1998) computed the first three radiative moment equations of Eq [15] and, by assuming that multipoles 
of /(n) higher than quadrupole vanish, found an otherwise exact form for the mean intensity perturbation, 5 J oc ik • 5/3, 
in a background with J — I, H = (e.g., such that 5 = A). Tacitly conforming with the premise of isotropy, wo re- write 
Eq [23] of Agol & Krolik (1998) by taking each of the perturbed diagonal components of the radiative field stress tensor, 
S K = (l/47r)(5 J di'dnnn7i/(n), to be identical and proportional to the perturbed radiative pressure Sprad- Thus SKu = ^SJ 
and, upon insertion of —5\nV in lieu of ik • df3, this yields the following relation for quasi-adiabatic perturbations of the 
radiative component of the fluid 

- --|(l-ia;) + — — }> — . (16) 



Prad 3l Ih U! 1 — iuj \ V 

Let us inspect the natural parameterization of this expression. When the non-difi^usive pressure is negligible, this index 

represents a truly adiabatic processes for wave modes only if Cj, h? &: fc^/iD <g; 1. Individually, both uj — a;(i!mfp/c) ~ 
|A|(c„/c)TjTjj^ — 2aA^(c„/c)^, and k ~ kT-Lr^^y. ~ 2y^|A|(c„/c), are expected to be small for the scales of interest in the 
magnetorotational instability problem. However, the ratio fe^/u; is not small and this has a non-trivial interpretation: 



k^ j_ . ^nilp/'''c"ddy i (c/Teddy) 



k2j w ^mfp |a;^,j,i l/c |A| 

teddy 



0{2i) (17) 



which identifles with the ratio of eddy diffusion rate to MRI growth rate. a4'fiv 
To zcroth order, from Eq [16] one thus has 

f^|(l + fP/<i)-i (18) 

which indicates that acoustic wave perturbations (i.e., with ui real) on the scale of the MRI are strongly damped if radiative 
pressure is the predominant source of stress. Indeed, Eq [18] may be interpreted as a form of energy equation characterizing 
radiative "heat conduction" out of compressive wave modes (Agol & Krolik 1998, compare this with Eq[51.12] of Mihalas & 
Mihalas 1984). 

If there is more parity between the diffusive and non-diffusive sources of stress, we flnd a generalized quasi-adiabatic 
index by defining a thermodynamic "quasi- volume" , V = , through Eq [16] as follows 

Prad 6 V 6 V 

while computing the diffusive stress contribution to thermodynamic processes in terms of the logarithmic differential of such 
quasi-volume: d\nV = i)d\nV. (This yields the effective volume filled by the radiative, "leaky" component of the fluid; the 
scale-dependent function fj is nothing but the logarithmic differential ratio between the effective and actual volumes occupied 
by the gas.) 

Wc consider two cases in turn: 1- an admixture of photon radiation plus an ideal gas (with pgas ~ Prad) and 2- an admixture 
of neutrinos (prad = 7/8ar*) plus gas/radiation/pairs (pgas = ll/12ar*) and where the latter species arc coUisionally coupled 
at all scales. Furthermore, for neutrino cooled flows we make the simpliflcation that the relatively cool neucleons (T s20 MeV) 
make a negligible contribution to the total pressure (see, e.g., Di Matteo et al. 2002) 

The standard lore (see, e.g., Mihalas & Mihalas 1984, Cliandrasekhar 1939) computation of the generalized adiabatic 
exponent, Fi, for an ideal gas plus a radiative component involves setting dQ = dU + dW = while working out expressions 
for dU & dW in terms of two logarithmic differentials, say, dlnT and dlnV. Using dQ = is clearly artiflcial but this is 
an adequate artifact when the main concern is with the non-elastic properties of the fluid and not with the amount of heat 
loss. Thus, we set dQtot = dQgas + dQrad = with both dQ gas ~ 0, and dQrad = (recall that the interactions between 
the two species are assumed to be entirely elastic). Computation of the specific heat contribution from the gas component 
is straightforward dQga,s/V = Pgas/(Fgas — l)dlnT -I- pgas dlnV; while use of the quasi-volume for the radiative component 
yields 

dQrad/t^ = 12praddlnT-|-4praddlnt> = 12praddlnT-h47jpraddlnV: (20) 

Elimination of dlnT in favor of dlnp is achieved by computing rfp(d In T, dlnV) for the combined gas. Furthermore, 
utilizing Fgas = 5/3 and the standard definition of Fi = —d[inv] Inp (= d[inp] Inp) yields 
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P _ (477/3 + l)(4/3 + l) + 12/3 + 3/2 
' (l+/3)(12/3 + 3/2) ' ^ > 

where /3 = Prad/Pgas- Naturally, the 5R(ri) agrees with the standard result (Chandrasekhar 1939, Mihalas & Mihalas 1984). 
Moreover, it is relatively straightforward to use the O'(ri) to find the decay rate of acoustic waves partially supported by 
radiative stress (aside from viscosity). 

For neutrino cooled accretion flows, Eq [20] is still valid under the understanding that Prad ~ 7/8aT^ corresponds to the 
neutrino pressure, while Pgas = ll/12aT* corresponds to photons and pairs and dQra.d/V = 12pgas dlnT + 4pgas din V. With 
these substitutions, an identical procedure to the one above yields 

_ 4 (l+r?/?) 



5 IMPACT OF RADIATIVE DIFFUSION AT MRI SCALES 

Thus far, wo have aimed to maintain the mathematical simplicity afforded by analytical solutions in order to gain physical 
insight on the phenomenology of radiative heat conduction and on its impact on the MRI. 

We have thus worked out solutions of a simplified version of the dispersion relation for non-axisymmetric MRI modes in 
the limit of purely horizontal fluid displacements. Since horizontal displacements maximize the efficiency of free-energy tapping 
from a difi^erential shear flow, this regime generally encompasses the unstable modes of fastest growth. The equation of motion, 
Eq [9], is applicable to standard heterogeneous fluids with arbitrary generalized index Fi. The ensuing dispersion relation, Eq 
[12], shows that the MRI is sensitive to a combination of two fluid properties-the "softness" of the fluid and the strength of 
the field-through a single parameter: A. Such parametrization, in turn, bears a direct connection to the compressibility of the 
modes through Eq [14]. Some physical insight may be gained by re- writing the compressibility parameter, 1 — A (oc Ap/p) as 
follows (c.f. Eq [14]) 

.A=|l+^^^'+^'" 



2 p^ 

while recapitulating that "2" represents the effective index of a magnetic field to perpendicular compression (see, e.g., Shu 
1992). MRI modes thus involve a large degree of compression when the effective sound speed of the material medium is much 
smaller than the Alfven speed. 

On the other hand, in §4 (c.f. Eq [17]) we have demonstrated that when the disk fluid is very hot such that a radiative, 
diffusive species (photons or neutrinos) constitutes a significant source of stress, radiative heat conduction precludes a formal 
treatment of the fluid as a single fluid at the scales of relevance to the MRI; that is, compression at these scales becomes 
non-adiabatic. A key conceptual point of this paper is that when the diffusive species dominates the stress tensor (e.g., 
radiation-pressure-dominated disks), the horizontal regime becomes inaccessible because of the brisk loss of pressure support 
from meridional diffusion^. Nevertheless, we are looking to read the gross properties of the fastest growth modes from the 2D 
dispersion relation Eq [12] by requiring parity among diffusive and non-diffusive sources of stress while finding an expression 
for the generalized quasi- adiabatic index of the combined fluid, Eq [21 & 22]. It is noteworthy that such parity among sources 
of stress is to be expected in neutrino cooled accretion flows such as those invoked in gamma-ray burst progenitor models. 

To place a zeroth order (recall [til ~ fe^ ^ 1, uj/k'^ ~ (^[Si]) upper limit on the quasi-adiabatic index, recall that for 
wave phase oc i(k • x — ujt) , the MRI growth rate corresponds to a purely imaginary wave frequency (i.e., since the wave phase 
is stationary at the co-rotation radius): w^^j^i = |. Thus, combining Eqs [17 & 18] yields 

^ = ^ 1 + ifeVK,,! ^ ^ + 1 {|} " ^) • 

An accurate estimate of the ratio (2^/|<iMnil arbitrary values of A(Fi, 0) may be found through the iterative use of 
Eq [13] (since 77 and thus Fi itself depends on g^/I^MRiDi but the value of the generalized index is rather insensitive to small 
changes in this ratio. For Keplerian rotation and in the weak field (i.e., incompressible) limit, we have 

4/I^MHil'^°5/4. 

On the other hand, the ratio k^/fcjj ':$> 1 for 2D motions so it should be clear that by setting > 2A;| we axe merely imposing 
a rough upper limit to the quasi-axiiabatic index: 

f<^W^°i (24) 



^ Note that the ordering of wavenumber components in the horizontal regime, k,f,kr<^kg, implies that the smallest length scale of the 
eddies is meridional; i.e., that the eddies are very nearly flat. 
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(again, the result we seek is rather insensitive to the precise value of ?7max)- 

Using ??max = 3/8, and prad = Pgas, in Eq [21], the generalized quasi-adiabatic index for an ideal gas plus photon radiation 
works out to be 

Fi — > 1+. (25) 
On the other hand, using Tjmax = 3/8, Prad = 7/8aT^ and pgas = ll/12aT"' in Eq [22], yields 

r, = ^ ^ 1- (26) 

1032 ' 

for a neutrtino cooled accretion flow. In cither case, the fluid behavior is close to isothermal (in a quasi-adiabatic sense) when 
the quEisi-adiabatic index associated with the diffusive species drops below the isothermal value F < 1. 



6 DISCUSSION 

Our phenomenological discussion is focused on two broad dynamical properties of the instability on the presumption that 
these help clarify the physical picture in the non-linear stage. Thus, the growth rate is identified in order of magnitude with the 
inverse of the correlation time of the turbulence and the geometrical regime of fastest growth, with the shape of the turbulent 
eddies. We also pay particular attention to the energy deposition process at the outer scale since a major portion of this 
paper is centered on the issue of entropy generation through non-adiabatic compression when a diflfusive species contributes 
non-trivially to the stress tensor. 

On the basis of the discussion in the previous section, we find that the root of a sluggish growth rate for the MRI 
(discovered by Blacs & Socrates 2001 for axisymmetric modes) when radiative pressure dominates the non-diffusive component 
is the adiabatic inaccessibility of the horizontal regime associated with fastest growth in standard fluids. For non-cixisymmetric 
modes, this regime becomes intrinsically non-adiabatic because of the brisk rate of heat transport out of anisotropic compressive 
perturbations. Because we begin with a toroidal field configuration where the non-axisymmetric modes are strongly compressive 
from the onset, the impact of radiative heat conduction is most dramatic and evidenced first hand from the dispersion relation. 

While in an ideal fluid the dominant eddies tend to have ke several times larger than kr and k^p, adding a radiative 
diffusive species will induce conductive losses that move mostly along Ig. This, in turn, means that the modes with high 
ke, i.e., horizontal modes, are preferentially damped. The magnitude of kg for the fastest growing modes is then be set by 
a marginal damping condition. Since the length scale for marginal damping is ~ fey (c.f. Eq 17), it follows that MRI eddies 
should be more nearly isotropic in such an environment, which goes along with the reduction in growth rate. For the same 
reeison, energy deposition on to the thermal bath (directly associated with the radiative species, and, secondarily, to the 
non-radiative species through scattering) occurs on the instability time scale. On the other hand, when the source of stress 
is equally divided among diffusive and non-diffusive particle species, the net effect is to soften the effective index of the fluid 
toward isothcrmality Fi 1 (c.f. Eq [25 & 26]) along with similar but more modest dynamical effects. 

These conclusions are broadly supported by the numerical simulations of Turner et al. (2001, 2003). Their report of 
a standard radiation-pressure-dominated a-disk, with prad S^Pgas — Pb shows that the non-linear outcome of the MRI is a 
porous medium with drastic density contrasts. Under nearly constant total pressure and temperature, the non-linear regime 
shows that density enhancements anti-correlate with azimuthal field domains-just as expected from the linear theory-and 
that turbulent eddies live for about a dynamical time scale while matter clumps are destroyed through collisions or by 
running through localized regions of shear on a similar time scale. Notably, the non-linear density contrasts may be quite 
large < pmax/pmin > X 10 when p^ad » Pgas- 

At a very fundamental level, clump formation is intimately connected to the effects of radiative heat conduction out of 
compressive perturbations. This can be understood by computing the polarization properties of the fastest growing modes 
in the horizontal regime: = — VA & j^^j <ti \£.r\, and reading from Eq [14] that at the linear stage of the instability 
there exists a converging fiow toward the Lagrangian displacement node of the modes (see Fig. 2 of Foglizzo & Tagger 1995). 
In a fluid with entirely elastic (adiabatic) properties, the pressure perturbation associated with such compression will act 
as a restoring force to de-compress the fluid in the non-linear stage. On the other hand, when the fluid has radiative heat 
conduction properties on the scale of the density perturbations, no such restoring force persists on time scales longer than the 
inverse of the growth rate (which, following Eq [17], is of the same order as the diffusion time). Material clumps thus formed 
will maintain their integrity and survive for longer times. When the diffusive and non-diffusive pressures are comparable, these 
arguments still hold but only in part. In this case, compressive perturbations in the non-linear stage will decompress only in 
proportion to the remaining pressure support from the non-diffusive species. 

Blaes & Socrates (2001) predicted qualitatively similar results for axisymmetric modes. Our analytical solutions agree 
with their result, |a)MRil ~ Cg/f^jj, in the limit kr = Cr ^ and Cg — > 0. This paper complements these findings and the 
non-linear numerical explorations of Turner et al. (2001, 2003) by presenting a simplified approximate dispersion relation 
for non-axisymmetric modes while interpreting the effects of radiative heat conduction in terms of an "ultra-soft" quasi- 
adiabatic index (i.e., F < 1) for the radiative contribution to the stress tensor. The quasi-adiabatic index and all derived 
thermodynamic quantities (including "quasi-specific" quantities that ensue from division by the quasi-volume, §4) are mere 
mathematical artifacts to manipulate thermodynamical relations in the usual manner. Furthermore, we have demonstrated 
that when there is parity among the sources of stress, the combined gas behaves isothermally. Our results are generally 
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applicable to optically thick, neutrino-cooled disks since the pressure contributions from photons, pairs and neutrinos, all 
have the same temperature dependence whereas only the neutrino component has radiative heat conduction properties on the 
time and length scales of the instability. 
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